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GRB Sky Distribution Puzzles 
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We analyze the randomness of the sky distribution of cosmic gamma-ray bursts. These events are 
associated with massive galaxies, spiral or elliptical, and therefore their positions should trace the 
large-scale structure, which, in turn, could show up in the sky distribution of fluctuations of the 
cosmic microwave background (CMB). We test this hypothesis by mosaic correlation mapping of 
the distributions of CMB peaks and burst positions, find the distribution of these two signals to 
be correlated, and interpret this correlation as a possible systematic effect. 



Key words: cosmology; cosmic microwave background 
observations — methods: data analysis 



gamma-ray bursts — cosmology: 



1. INTRODUCTION 

The publication of comparatively large catalogs of 
gamma-ray bursts based on the results of the surveys 
performed within the framework of space missions 
of the Italo-Dutch BeppoSAX satellitq^ (Satellite per 
Astronomia X, and Beppo in the honor of Giuseppe 
Occhialini) [1], and NASA's Compton observatory, 
where the BATSE experimenJl (Burst and Transient 
Source Experiment) [2] was performed, made it pos- 
sible to study the properties of the space distribution 
of gamma-ray bursts, which, to a first approximation, 
are uniformly distributed in the sky. 

Among recent studies concerning this subject we 
distinguish the analysis of the possible relation be- 
tween the distribution of gamma-ray bursts from 
the BATSE catalog and the large-scale structure by 
Williams and Frey [3] . The authors analyze the effect 
of the local large-scale structure on the apparent posi- 
tions of gamma-ray bursts in the case of weak lensing. 
This effect may show up for distant (z > 4) events 
as an anticorrelation of the positions of gamma-ray 
bursts and galaxies in clusters. The authors studying 
gamma-ray bursts report such an anticorrelation for 
galaxies at z ~ 0.2 — 0.3, which they found based 
on the optical magnitudes and positions of galaxies 
measured using the APM (Automatic Plate Measur- 
ing Facility). 

In another paper, Meszaros et al. [4] use vari- 
ous methods (the Voronoi tesselation diagrams, the 



^ http : //www. asdc . asi . it/bepposax/ 

^ jhttp : //www.batse .msf c .nasa. gov/batse/| 



minimal spanning tree, and multifractal spectra) to 
analyze the distribution of the BATSE gamma-ray 
bursts. The above authors break the list of objects 
into subsamples of sources with different signal dura- 
tions {t < 2 s, 2 < i < 10 s, t > 10 s) and apply the 
above methods to each subsample. They found the 
results for the first two groups to deviate from the 
results obtained for uniformly distributed simulated 
data. Based on the data obtained the authors [4 dis- 
cuss the validity of the cosmological principle. 

We use the BeppoSAX (781 source, the energy 
interval covered: 0.1-200 keV) and BATSE (2037 
sources, 20keV-2MeV) catalogs adopted from the 
sites of the corresponding experiments. We subdivide 
each catalog into two subsamples containing short 
(i.e., lasting t < 2s) and long {t > 2s) events. We 
thus obtain four source subsamples, which we ana- 
lyze separately using the same procedure. Figure 1 
shows the positions of all gamma-ray bursts from the 
BeppoSAX and BATSE catalogs. Figure 2 shows the 
positions of short and long BeppoSAX and BATSE 
gamma-ray bursts on separate sky charts. 

In this paper we analyze the statistical correla- 
tion properties of the GRB distribution on the sky 
with respect to the distribution of the cosmic mi- 
crowave background (CMB). We proceed from the 
assumption that gamma-ray bursts should be asso- 
ciated with massive galaxies, and that the positions 
of these galaxies are correlated with the large-scale 
structure. In this case there is hope that we may find 
deviations in the sky distribution of CMB fluctuations 
in the areas where the bursts are projected. Note that 
the large-scale structure may show up in CMB flue- 
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Figure 1: Distribution of gamma-ray bursts on the sky. The gray and white symbols (the circles and crosses 
correspond to short and long bursts, respectively) show the BeppoSAX and BATSE data, respectively. 



Figure 2: Distribution of the gamma-ray bursts on the sky. The top left panel shows the BeppoSAX data for 
bursts with t < 2 s. The top right panel shows the BeppoSAX data for bursts with t > 2 s. The bottom left panel 
shows the BATSE data for t < 2 s. The bottom right panel shows the BATSE data for t > 2 s. 



tuations via the Zel'dovich-Sunyaev effect and the 
integrated Sachs- Wolfe effect [B]. This former effect 
manifests itself via the interaction of CMB photons 
with hot electrons of galaxy clusters and is ob- 
served on the angular scales of about 4' or smaller. 
The latter effects occurs when the photons move in a 
gravitational field with a variable potential that arises 
during the formation of the large-scale structure and 
expansion of the Universe, and should be observed on 
the angular scales greater than 200'. 

To study the correlation properties of the posi- 
tions of gamma-ray bursts and CMB peaks, we use 
the CMB map based on the results of the experi- 
ment that has been conducted for seven years within 
the framework of the WMApjl (Wilkinson Microwave 
Anisotropy Probe) experiment aimed to integrate the 
signal from the entire sky [7]. The CMB signal was 
reconstructed from multifrequency observations us- 



^ |http : //lambda.gsf c .nasa. gov] 



ing the method of internal linear combination (ILC) 
of the background components [S]. This procedure 
yielded the CMB map, which is also referred to as 
the ILC map, and which is used to analyze low-order 
harmonics with multipole numbers £ < 150. The ILC 
map is based on observations made in five channels: 
23 GHz (the K band), 33 GHz (the Ka band), 41 GHz 
(the Q band), 61 GHz (the V band), and 94 GHz (the 
W band). 

Given that one of the central problems in statis- 
tical studies of gamma-ray bursts is due to the large 
error box size, about l°xl°, of the available source 
positions, we operate with pixels of about the same 
or even greater size in order to avoid uncertainties in 
the analysis of the sky distribution of the events. 

The methods we use in this paper and our results 
are laid out in the following way. First we describe 
the methods that we use to pixelize the distribution of 
bursts and to correlate the maps. We then analyze the 
statistics of CMB deviations in the gamma-ray burst 
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areas. At the next stage, we compare the BATSE- 
BeppoSAX and BATSE-CMB correlation maps. We 
finahy discuss the resuhs in the conchiding section. 



2. MOSAIC 
METHOD 



CORRELATION 



To analyze the map properties on different angular 
scales, we expand the signal distributed on the sphere 
into the spherical harmonics (multipoles) : 
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Here AS{6, (j)) are the variations of signal on the 
sphere in polar coordinates; £ is the number of the 
multipole, and m is the number of the mode of the 
multipole. The spherical harmonics are determined as 
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where Pp{x) are the associated Legendre polynomi- 
als. The expansion coefficients aim for a continuous 
function AS{x, </>) can be written as: 



air, 
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dx / AS{x,(j))Yl^{x,(t))d(t) , 
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where Yg^^ denotes the complex conjugate of Yim- 

The correlation properties of two maps of the sig- 
nal distribution on the sphere can be described, on 
a given angular scale, by a correlation coefficient for 
the corresponding multipole i as: 



K{i) = - 



(4) 



where tem and sim are the variations of two signals in 
a harmonic representation. The coefficient K{£) can 
be used to assess the correlation between the harmon- 
ics on the sphere, i.e., to compare the properties of 
maps on a given angular scale. However, in the case 
of a search for correlated areas, which do not repeat 
in other regions of the sphere, this approach smears 
such single areas in the process of averaging over the 
sphere within a certain harmonic. In this case it be- 
comes practically impossible to identify the correlated 
areas. 

Verkhodanov et al. [3] proposed an approach, 
which was implemented in the second release of the 
GLESP package [T^ (the difmap utility). The pro- 
posed procedure makes it possible to find correla- 
tions between two maps in the areas of a certain 



angular size. In this method, each pixel with num- 
ber p subtending the solid angle Sp is assigned the 
cross-correlation coefficient between the data of the 
two maps on the corresponding area. Thus a cor- 
relation map is constructed for two signals T and 
S, where the value of each pixel p [p — l,2,...,7Vo) 
and Nq is the total number of pixels on the sphere) 
with the subtending angle 5p, and computed for the 
sphere maps with the initial resolution determined by 
(■max is equal to 



E E iTi9„<j>,)-TiEp))iSie„^,)-S{Ep)) 



(5) 



Here T{6i, ipj) is the value of the signal T in the pixel 
with the coordinates {9i,(f)j) for the initial resolution 
of the pixelization of the sphere; S{9i, (f>j) is the value 
of the other signal in the same area; T(Sp) and >5'(Sp) 
are the mean values averaged over the area Sp and 
obtained from the data of the maps with higher res- 
olution determined by £max , and and as^ are the 
corresponding standard deviations in the area consid- 
ered. 



3. CMB SIGNAL STATISTICS IN 
THE AREAS OF GAMMA-RAY 
BURSTS 

In the first method the statistical properties of dis- 
tributions of GRB positions versus the CMB signal 
are tested by simply analyzing the pixel value mea- 
surements on the CMB maps. We performed pixel 
measurements using the mapcut procedure of the 
GLESP package. Figures 3 and 4 show the positions 
of gamma-ray bursts of the BeppoSAX and BATSE 
catalog subsamples on the CMB maps with the res- 
olution of 260' {(max = 20) and 36' {£max = 150). 
We chose the resolution of these maps to match the 
expected angular scale of the Sachs- Wolfe effect and 
the limiting resolution of the WMAP mission. 

To search for the eventual correlations, we com- 
puted the number of GRB positions that fall onto the 
CMB pixels with negative signal fluctuations (which 
may be due to the effects described above) in the 
CMB maps of different resolution. The Table lists the 
statistics of CMB pixel values in the gamma-ray burst 
areas for the long- and short-event BATSE and Bep- 
poSAX subsamples. It includes the total number of 
sources in subsamples; the number of sources located 
in CMB pixels with negative fluctuation values; the 
expected number of pixels with negative CMB ampli- 
tudes according to the results of 200 random Gaus- 
sian CMB signal simulations modeled in terms of the 
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Figure 3; Positions of gamma-ray bursts from different suhsamples on the CMB maps with a resolution of 
imax — 20. The top left and top right panels show the BeppoSAX data for bursts with t < 2 s and t > 2 s, 
respectively. The bottom left and bottom right panels show the data for the BATSE bursts with t < 2s and 
t > 2 s, respectively. 




Figure 4: Positions of gamma-ray bursts from different .suhsamples on the CMB maps with a resolution of 
(■max = 150. The top left and top right panels show the BeppoSAX data for bursts with t < 2 s and t > 2 s, 
respectively. The bottom left and bottom right panels show the data for the BATSE bursts with t < 2s and 
t > 2 s, respectively. 



ACDM cosmology, and the Icr-scatter of these quan- 
tities. 

Figures 5 and 6 show the distributions of the CMB 
fluctuations for four gamma-ray burst suhsamples 
and two CMB maps with different resolutions. The 
dashed lines show the expected Icr scatter of CMB 



values in the ACDM cosmological model. The bot- 
tom left panel in Fig. 6, which shows the distribution 
of fluctuations on the CMB map with a resolution of 
(■max — 150 in the areas of short gamma-ray bursts 
of the BATSE catalog, stands out conspicuously. In 
the positive part of the plot there is a peak, which 
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Figure 5: Distribution of CMB fluctuations in the WMAP pixels corresponding to the GRB positions on the 
maps with a resolution of imax = 20. The top left and top right panels show the BeppoSAX data for the bursts 
with t < 2s and t > 2s, respectively. The bottom left and bottom right panels show the data for the BATSE 
bursts with t < 2 s and t > 2 s, respectively. The dashed lines show the la scatter of CMB values in the A CDM 
cosmological model. 



makes the distribution deviate from Gaussianity. Fig- 
ure 7 shows the positions of these bursts superim- 
posed onto the CMB quadrupole. The statistical sig- 
nificance of this feature, which we estimated by gen- 
erating 10000 simulated CMB realizations with the 
power spectrum corresponding to the ACDM model, 
is equal to 7 x 10^**. 

To analyze the distribution of GRB positions cor- 
responding to the excess peak, we used the GLESP 
software package jlO| to pixelize the map of burst 
positions. We chose the pixel size 700' x 700' in or- 
der to make the maximum pixel value (the number 
of events inside the corresponding area) no less than 
three, and to provide a significant dynamic range for 
the harmonic analysis. Figure 8 shows the map of 
selected bursts pixelized in such a way. The distribu- 



tion of events on the sphere shown in the image is 
clearly non-uniform: it is concentrated near the eclip- 
tic and/or equatorial poles, which is immediately ap- 
parent in the smoothed map (Fig. 9). Figure 9 also 
shows the ecliptic and equatorial coordinate grids, 
and it is evident from this figure that the asymme- 
try shows up both in the distribution of the signal 
power with respect to the equatorial plane, and in 
the number of spots that concentrate in the South- 
ern Hemisphere in both coordinate systems. Note also 
that the hot spot above the Galactic center lies in the 
equatorial plane. 

To reveal the asymmetries more clearly, we 
demonstrate a quadrupole of the smoothed map 
(Fig. 10) with the ecliptic and equatorial coordinate 
grids superimposed, as we do it in Fig. 9. It is evident 
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Figure 6: Distribution of the CMB fluctuations in the WMAP pixels corresponding to GRB positions on the 
maps with a resolution of i^nax = 150. The top left and top right panels show the BeppoSAX data for the bursts 
with t < 2s and t > 2s, respectively. The bottom left and bottom right panels show the data for the BATSE 
bursts with t < 2 s and t > 2 s, respectively. The dashed lines show the la scatter of CMB values in the A CDM 
cosmological model. 




Figure 7: Positions of short BATSE bursts corresponding to the excess in the histogram in Fig. 6 (the bottom 
lefl panel) on the map of the quadrupole ILC component. 
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Figure 8: Pixelized map of positions of short gamma-ray bursts of the BATSE catalog corresponding to the 
excess in the histogram in Fig. 6 (the bottom left panel). The pixel size lOCf x 70C/ is chosen in a way to make 
the maximum pixel value — the number of events inside the corresponding sky area — greater than or equal to 
three. 



Figure 9: Smoothed sky map corresponding to the pixelization used in Fig. 8, with the ecliptic (right) and 
equatorial (left) coordinate grids superimposed. 



Table 1: Table. Statistics of the CMB pixel values in 
the GRB areas for the BATSE and BeppoSAX sub- 
samples. The columns are: t, the duration (s); mis- 
sion name; resolution of the CMB map ( the number 
of the multipole); total number of gamma-ray sources 
in the subsample (Nt); the number of sources (N^.) 
located in the CMB pixels with negative fluctuation 
values; expected number of pixels with negative CMB 
amplitudes according to the results of 200 random 
Gaussian CMB signal simulations, modelled in terms 
of the ACDM cosmology, and the la -scatter of these 
quantities 



t, s 


Mission 


^max 


Nt 




Model 


<2 


BATSE 


150 


497 


244 


249±11 


>2 


BATSE 


150 


1540 


763 


769±19 


<2 


BATSE 


20 


497 


250 


248±13 


>2 


BATSE 


20 


1540 


781 


768±32 


<2 


BeppoSAX 


150 


87 


43 


44±5 


>2 


BeppoSAX 


150 


694 


339 


347±15 


<2 


BeppoSAX 


20 


87 


50 


44±5 


>2 


BeppoSAX 


20 


694 


346 


348±18 



from the positions of quadrupole spots that: (1) the 
poles of both coordinate systems are located on the 
outskirts of hot spots; (2) cold spots are located sym- 



metrically with respect to the equatorial planes, and 
(3) the Galactic center is located in the saddle-shaped 
area between cold and hot spots. 

Our analysis of gamma-ray burst positions in the 
area of the peak in the distribution of short BATSE 
bursts with respect to the CMB revealed their un- 
expected sensitivity to local (near-Earth) coordinate 
systems. To analyze this problem in more detail, we 
pixelize all the four subsamples of the GRB catalogs 
using the method of mosaic correlation described in 
the previous section. 

4. GRB AND WMAP ILC DISTRIBU- 
TION CORRELATION MAPS 

To verify and refine the correlation properties of the 
maps of GRB positions and CMB fluctuations, we 
pixelized the maps of GRB positions for four gamma- 
ray burst subsamples (Fig. 11). Like at the previous 
stage, we chose the pixel size 200' x 200' {imax = 26) 
in order to ensure that the maximum pixel value — the 
number of events in the corresponding area — would 
be not less than three. 

We then performed the mosaic correlation for 
the BATSE-BeppoSAX and BATSE-CMB pairs 
(Fig. 12). The pixel size for these correlations was 
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Figure 10: A quadrupole of the smoothed map (Fig. 9) with the ecliptic (right) and equatorial (left) coordinate 
grids superimposed. 




Figure 11: Pixelized maps of event positions in gamma-ray hurst subsamples. The pixel size is 20ff x 200'. The 
top left and top right panels show the BeppoSAX data for bursts with t < 2 s and t > 2 s, respectively. The 
bottom left and bottom right panels show the data for the BATSE bursts with t < 2 s and t > 2 s, respectively. 



equal to 500' x 500'. 

We can now use the expansion defined by for- 
mula ([1]) to compute the angular power spectrum 
of the map: 



C{£) 



1 



2i + i 



\aio\ 



2E 



(6) 



The power spectrum allows us to identify the main 
harmonics contributing to the correlation map. Fig- 
ure 13 shows the power spectra of the correlation coef- 
ficient maps (mosaic correlation) between the BATSE 
and CMB data {imax = 26) and between the BATSE 
and BeppoSAX gamma-ray positions. 

Three correlation power spectra stand out among 
other correlations: those between the positions of 
short BATSE and BeppoSAX bursts (the top left 
panel in Fig. 13) and between short and long BATSE 



bursts and CMB fluctuations (top and bottom right 
panels, respectively in Fig. 13). The peaks of these 
spectra exhibit features that formally lie beyond the 
2a level, and, in the case of the harmonics i = 1, 
4, 7, 9, even beyond 3(t. To estimate the statistical 
significance of these deviations, we performed numer- 
ical simulations of the zero hypothesis by generating: 

(1) 10,000 random realizations of the set of sky posi- 
tions of the BeppoSAX gamma-ray bursts assuming 
that they are uniformly distributed in the sky and 

(2) random CMB maps with the Gaussian amplitude 
distribution corresponding to the ACDM cosmology, 
and computed the correlations and the correspond- 
ing power spectra for these simulated data in the 
same way as we did it for real data. Figure 14 shows 
the probabilities for a simulated realization to reach 
the deviation levels of the features observed in Fig. 13 
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Figure 12: Correlation maps for the CMB data (imax — 26) and gamma-ray burst positions in the Galactic 
coordinate system. The pixel size for these correlations is equal to 50Cf x 500'. The top left panel shows the 
correlation between BATSE and BeppoSAX data, t < 2 s, with the ecliptic coordinate grid superimposed. The 
top right panel shows the correlation between BATSE (t < 2 s) and CMB data. The bottom left panel shows the 
correlation between BATSE and BeppoSAX data, t > 2 s. The bottom right panel shows the correlation between 
BATSE (t>2s) and CMB data. 



within the framework of the zero hypothesis i.e., the 
significance levels of the peaks in the power spec- 
trum. The resulting significance levels for the mul- 
tipoles £ = 1, 4, 7, and 9 of the correlation maps 
between the short BATSE and BeppoSAX bursts are 
equal to 0.0011, 0.0012, 0.0001, and less than 0.0001, 
respectively. At the same time, the significance level 
for the octupole of the correlations between the long 
BATSE and BeppoSAX bursts is equal to 0.0325, and 
that for the quadrupole of the correlation map be- 
tween the long BATSE bursts and CMB fluctuations 
is equal to 0.0614. Figure 15 shows the map of one 
of the harmonics (£ = 7) in the distribution of the 
correlation coefficients between the positions of short 
BATSE and BeppoSAX bursts. Figure 16 shows the 
maps of the peak harmonics (£ = 3 for short and 
£ = 2 for long bursts) in the spectrum of the maps of 
correlations between the BATSE and CMB. 

A distinguishing feature of the maps shown with 
superimposed coordinate grids are the positions of 
the poles of coordinate systems. It is evident from 
Fig. 15 that the ecliptic poles are located at spe- 
cial points — in the saddles between the maxima and 
minima of the signal distribution on the map. Sim- 
ilarly, the equatorial poles are located in the sad- 
dles in the map of the correlation octupole (the left 
panel in Fig. 16). We can see in the right panel (the 
quadrupole) of Fig. 16 that the poles are located in 



the minima of the quadrupole. We use the method 
described above based on the generation of 10,000 
random realizations of the GLESP pixelization with 
102 pixels at the equator to estimate the statistical 
significance of such a configuration i.e., the proba- 
bility for the minima of the quadrupole to occur in 
the 5-degree radius areas centered on the equatorial 
poles, and find it to be equal to 0.0035. 

5. DISCUSSION 

Our analysis of correlations between the maps of GRB 
positions and the CMB shows that such correlations 
do exist. Of great interest is the orientation (phase 
properties) of correlation maps, where we can point 
out two points: (1) in the case of the correlations be- 
tween the short GRB positions in different catalogs 
and their correlations with the CMB, the features 
found (the positions of the poles) are observed both in 
the equatorial and ecliptic coordinate systems; more- 
over, despite a small number of short BeppoSAX and 
BATSE events, the correlations between their posi- 
tions in 500' X 500' windows exhibit a chain of events 
in the ecliptic plane (the top left panel in Fig. 12), 
and a predominant occurrence of the correlated pixels 
in the Southern Hemisphere; (2) the correlations be- 
tween the long BATSE events and CMB fluctuations 
exhibit features in the equatorial coordinate system, 
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Figure 13; Power spectra of the correlation coefficient maps between the BATSE GRB positions and BeppoSAX 
and CMB data (£max = 26j (the solid line). The pixel size of the correlations is SOff x500' . The top left panel 
shows the spectrum of correlation between the BATSE and BeppoSAX data, t < 2 s. The top right panel shows 
the spectrum of correlation between the BATSE (t < 2 s) and CMB data. The bottom left panel shows the 
spectrum of correlation between the BATSE and BeppoSAX data, t > 2 s. The bottom right panel .shows the 
.spectrum of correlation between the BATSE (t > 2 s) and CMB data. The dashed lines show the ±cr, ±2a, and 
±3(7 deviations from the mean, where a is estimated from the .spectra of 200 simulated random realizations of 
the zero hypothesis, which assumes a uniform distribution of gamma-ray bursts in the sky and the Caussian 
CMB amplitude distribution corresponding to the ACDM cosmology. 



while the joint probabihty of high quadrupole amph- 
tude and of the quadrupole minimum occurrence at 
the equatorial pole is practically equal to zero in the 
case of correlation with random maps. 

The correlations found between the GRB posi- 
tions and the CMB, that are sensitive to the equa- 
torial coordinate system must be due to the system- 
atic effects. This way, a relatively greater number 
of gamma-ray events near the equatorial poles may 
be due to longer exposures of the satellite cameras 
due to the observing method employed. Both satel- 
lites that observed the gamma-ray bursts moved in 
rather low Earth orbits, rendering the areas located 
near the celestial equator periodically unobservable. 
At the same time, for the CMB data such a sensitiv- 
ity to the equatorial coordinate system is impossible 



to explain in terms of such a simple model, as the 
microwave background data were obtained onboard 
the WMAP satellite, which rotates about the L2 La- 
grangian point. Note that in our previous paper [11] 
we also found some correlations to be "aware" of the 
equatorial system. We do not rule out the possible 
contribution from the Earth's magnetic field, which 
has a large extent and shows up as large-scale corre- 
lations of the microwave background. However, we do 
not understand the mechanism of such correlations. 

Neither do we understand the possible interrela- 
tion between the ecliptic plane and the positions of 
gamma-ray bursts. Here the situation is reversed: the 
ecliptic features in the CMB data have already been 
discussed recently [12l [131 [11] , however, further stud- 
ies are needed to understand what may link gamma- 



248 



O.V. Verkhodanov,V.V. Sokolov, M.L. Khabibullina, S.V. Karpov 




Figure 16: Maps of the harmonics identified in the power spectrum (Fig. 13). Correlation maps between the 
BATSE hurst positions (t > 2 s) and the BeppoSAX data (the top panel, £ = 3) and between the BATSE burst 
positions (t > 2 s) and the CMB (the bottom panel, 1 = 2) with the equatorial coordinate grid superimposed. 




2 4 6 a 10 

Multipole L 



Figure 14: Significance levels for the deviations from 
the mean level for the power spectra of the correlation 
maps: short gamma-ray burst BATSE-BeppoSAX 
correlations (the thin dashed line); long gamma- 
ray burst BATSE-BeppoSAX correlations (the thick 
dashed line); short BATSE bursts and CMB fluctu- 
ations (the thin solid line), and long BATSE bursts 
and CMB fluctuations (the thick solid line). The sta- 
tistical significance was estimated via direct numer- 
ical simulations involving 10000 random realizations 
of the zero hypothesis, which consists in the uniform 
GRB distribution in the sky and the Gaussian CMB 
amplitude distribution corresponding to the ACDM 
cosmology. 




Figure 15: Map of the seventh harmonic (£ = 7) in 
the expansion of the map of correlations between the 
positions of short BATSE and BeppoSAX bursts with 
the equatorial coordinate grid superimposed. 



ray evens and the ecliptic plane. One of the possible 
hypotheses explains this feature as a selection effect 
due to the fact that the observing instruments of the 
satellites are turned away from the Sun, and this may 
put the plane of the ecliptic in a special position. 

Note also that the correlation properties of the 
CMB, which show up in the "near-Earth" coordinate 
systems when the data with "randomly" distributed 
events are used, are indicative of the non-Gaussian 
nature of the data for low-order multipoles, which 
may be due either to the systematic effects, or to a 
hitherto unexplored effect in the near-Earth space. 
We consider further studies of this correlation to be 
of especially great interest given that the new high- 
quality data are expected to be provided by the Fermi 
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and Planck missions. 
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